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Abstract 

Background The recent Jefferson Lab data for the nuclear transparency in electroproduction 
have the potential to settle the scale for the onset of color transparency (CT) in vector meson 
production. 

Purpose To compare the data to calculations in a relativistic and quantum-mechanical Glauber 
model and to investigate whether they are in accordance with results including color trans- 
parency given that the computation of p-nucleus attenuations is subject to some uncertain- 
ties. 

Method We compute the nuclear transparencies in a multiple-scattering Glauber model and ac- 
count for effects stemming from color transparency, from p-meson decay, and from short- 
range correlations (SRC) in the final-state interactions (FSI). 

Results The robustness of the model is tested by comparing the mass dependence and the hard- 
scale dependence of the A{e, e'p) nuclear transparencies with the data. The hard-scale de- 
pendence of the (e, e'p^) nuclear transparencies for ^^C and ^^Fe are only moderately affected 
by SRC and by p^-decay. 

Conclusions The RMSGA calculations confirm the onset of CT at four-momentum transfers of a 
few (GeV/c)^ in p meson electroproduction data. A more precise determination of the scale 
for the onset of CT is hampered by the lack of precise input in the FSI and /o-meson decay 
calculations. 
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I. INTRODUCTION 



Color transparency (CT) is a color coherence effect that emerges from perturbative quan- 
tum chromodynamics in exclusive knockout reactions at high four-momentum transfers 
[H 12]. In these reactions, a hadron is produced in a small-sized configuration (SSC) with 
all quarks sitting close together in the transverse plane. The color interactions with the 
surrounding nuclear medium cancel each other and the hadron can propagate unattenuated 
as the common final-state interactions (FSI) between the tagged hadron and the nuclear 
environment vanish. The SSC can also be produced in nonperturbative conditions [SlH]. In 
this regime, the SSC will evolve to its stable hadronic state over a certain formation length 
//. During its formation the hadron is subjected to reduced interactions with the nuclear 
medium. In order to observe a CT effect under those conditions, the formation length should 
be of the order of the nuclear radius (// ~ Ra)- Observation of the onset of CT at a cer- 
tain energy scale can teach us about the cross-over point between ordinary nuclear matter 
and quark-gluon matter. The identification of this transition point is of great importance 
to nucleon-structure studies, as CT is a necessary condition for the validity of the QCD 
factorization theorems which are commonly applied when interpreting the data |5l E] . 

The measured observable in search for CT is the nuclear transparency T, defined as the 
ratio of the cross section per target nucleon for a process on a nucleus to the cross section 
of the process on a free nucleon (T = a^lA(j^\ Accordingly, the nuclear transparency 
provides a measure of the integrated attenuation of the nuclear medium on the tagged 
hadrons in some (semi-)exclusive reaction. One can study the hard-scale dependence of 
the transparency for a certain target nucleus A, or the A dependence at a fixed value of 
the hard-scale parameter. If CT effects were to appear at a certain value of the hard-scale 
parameter, the nuclear transparency would be observed to overshoot the predictions from 
traditional nuclear-physics calculations. The measurement of the onset and magnitude of 
the CT effect allows one to constrain models describing the evolution of a SSC into a hadron. 

Experimentally, CT effects have been observed in the measurement of the cross section 
of diffractive dissociation of 500 GeV/c pions into dijets in the E791 experiment at Fermilab 
[7]. At intermediate energies, no sign of CT was observed in A{e,e'p) measurements on 
a variety of nuclear targets and four-momentum transfers ^ 8 (GeV/c)^ pHT5J- The 
nuclear ^^C(p, 2p) transparencies were studied at Brookhaven National Laboratory (BNL) 
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[T5HT7] . The transparency first shows a rise with increasing incoming proton momentum 
and drops for momenta larger than 9 GeV/c. This proton-momentum dependence is at 
odds with traditional nuclear-physics calculations predicting ^^C(p, 2p) nuclear transparen- 
cies which are more or less constant with proton momentum. The BNL ^^C(p, 2p) results 
are currently not considered as a clean sign of CT, and competing effects stemming from 
nuclear filtering [19] or from threshold mechanisms for charm resonance production [20] 
have been proposed to explain the observations. 

In recent years, several experiments have measured the transparencies in semi-exclusive 
meson production reactions. As a meson is a more compact object than a baryon, it should 
be more likely to produce a meson SSC and observe the onset of CT at intermediate en- 
ergies. Two recent Jefferson Lab (JLab) experiments that measured the transparency of 
pions in photo- and electroproduction [211 E2] consistently agreed with various independent 
calculations provided that CT effects are included [231426] . More recently, a JLab Hall-B 
experiment has measured nuclear transparencies in semi-exclusive electroproduction on 
^^Fe and targets [27]. A gain, the data agree favorably with calculations including CT 
effects [281 EH]. These results strengthen the case for an onset of CT at four-momentum 
transfers of a few (GeV/c)^ in exclusive meson production reactions. For a recent review of 
the CT phenomenon, see Ref. [30] . 

There are issues in nuclear p transparency calculations, like the absence of pN scattering 
data, which complicate the interpretation of the calculations and induce some uncertainties. 
It is one of the purposes of this paper to investigate theses issues in more detail and to 
study the robustness of the computed nuclear p transparencies. We compare the recent 
A{e, e'p^) transparency data with calculations in a relativistic multiple-scattering Glauber 
approximation (RMSGA) model. The RMSGA model has been used successfully to predict 
nuclear transparencies for A{e,e'p) reactions [311 [32], for A{p,2p) reactions [331 El], for 
^4(7, TT+p) and A (e, e'7r+) reactions [231 [21] and for quasi-elastic neutrino-induced processes 
[35] . In Sec. |ll] we sketch the RMSGA model while highlighting some important issues 



emerging for the p° nuclear transparency calculations. Numerical results are shown in Sec. Ill 



for kinematics of the JLab experiment. Sec. IV states our conclusions. 
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II. MODEL 



a. Relativistic multiple-scattering Glauber approximation Finding its roots in optics, 
Glauber multiple-scattering theory [35] describes the small-angle scattering of particles in 
the eikonal approximation. Thereby the scattered wave function is a plane wave multiplied 
with an eikonal phase. The eikonal approximation is valid when the wavelength of the par- 
ticle is small in comparison with the typical interaction range of the scattering particles. 
For baryons and mesons this criterion translates into momenta higher than a few hundred 
GeV/c. In Ref. [31] we have introduced a relativistic version of Glauber multiple-scattering 
theory. As the helicity conserving amplitude is assumed to dominate hadron-nucleon high 
energy scattering, the Glauber eikonal phase is a scalar in spin space. To describe multiple- 
rescattering the frozen approximation is used and the individual eikonal phases are multi- 
plied. In the RMSGA model, the eikonal Glauber phase at some spatial point {b, z) reads 
for the ejected p° 

QiP,z)= n j dr'\(t)^{r')\^[l-e{z-z')V,N{b' -b)\ . (1) 

Here, the (6, z) coordinate system has its z-axis along the momentum, 0Q,(r ) are the Dirac 
single-particle wave functions of the residual nucleons with quantum numbers a = nurrij 
obtained from the Serot-Walecka model [3Z] , and ctj characterizes the nucleon on which the 
p° is created. The profile function F^at in Eq. M is commonly described by a Gaussian 



where a*'^, epN, PpN are energy dependent parameters connected to pA^ scattering. Due to 
the lack of data for pA^ scattering it has become customary to use educated estimates for the 
CpjV) ^pN^ PpN [38] based on the corresponding values for vrA^ scattering. In vrA^ scattering the 
parameters display very little energy dependence for pion lab momenta p-,^ > 1.5 — 2 GeV/c 
and it has become customary to adopt energy-independent parameters for pA^. In this work, 
we study the impact of these uncertainties for the computed nuclear p transparencies. 

b. Short-range correlations The RMSGA model accommodates the possibility to in- 
clude short-range correlations (SRC) in the modeling of the FSI. We use the information 
that a nucleon is present at the spatial point of the hard interaction. Due to its finite size, 
the presence of the nucleon induces local fluctuations in the nuclear density. The inclusion 



of SRC in the FSI is technically achieved in the following way [21]. First, the squared single- 
particle wave functions in Eq. ([T]) can be connected to the one-body density of the target 
nucleus pj|^(r) (normalized as / dfp^^{r) = A) 

Here, is the ground-state wave function of the target nucleus, obtained by anti- 

symmetrizing the product of the single-particle wave functions (pa- Even in a small nu- 
cleus like '^He, the above approximation (better known as the thickness approximation) 
marginally affects the predicted effect of FSI [31]. In a second step, the averaged den- 
sity p^^ (r ) can be substituted with the ratio of the two-body density p^^ (normalized as 
/ dfi J df2p^'^{fi,f2) = A{A — 1)) and the one-body density: 

where n is the spatial coordinate corresponding with the hard interaction. One can include 
SRC in the two-body density by adopting the functional form |39j : 

p[J'(fi,f2) ^ ^^^1 (r i) P A (r i) p^A^r 2)1 (r 2) g{r 12) , (5) 

where (7(^12) is the Jastrow correlation function [10] and 7(r ) a function that imposes the 
normalization of the two-body density obtained as the solution of an integral equation. With 
the above expression for the SRC-corrected two-body density, Eq. ^ becomes 

p5'(f) ^ 7(f )p5^(f)7(fi)^(|f -fi|) =pf{f,fi). (6) 

In summary, the FSI calculations can be corrected for SRC by replacing |0Q,(r )|^ with 
p^(r, ri)/^ in Eq. . To illustrate the effectiveness and robustness of the RMSGA model 
in describing a variety of knockout reactions, we show in Fig. [T] RMSGA predictions for 
the A{e, e'p) nuclear transparencies and compare them to the published data. Neither the 
calculations nor the data include the commonly applied (and rather arbitrary) correction 
factors ca that account for shifts in the nuclear spectral function to higher missing momenta 
and energies due to SRC. For a discussion on this subject, we refer to Ref. [30] • We observe 
an excellent agreement between the RMSGA calculations and the world data on A{e, e'p) 
nuclear transparencies. Both the hard-scale dependence and the mass dependence of 
the data are nicely predicted. It is noteworthy to mention that the A (e, e'p°) calculations 
presented below use identical nuclear-structure input. 
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FIG. 1: (Color online) Nuclear transparencies versus for A{e,e'p) reactions in quasi-elastic 
kinematics. Solid black lines show the RMSGA calculations corrected for SRC according to the 
replacement of Eq. Q in Eq. ([T]). Data are from Refs. f8] (black squares), O HO] (red hexagons), 
|12j (green stars), [H] (cyan diamonds) and I13j (blue circles). Data and calculations do not 
include the ca factor applied in [32] . 



Color transparency The effects of color transparency are implemented by means of 
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the quantum diffusion model of Ref . [H] . Thereby, the position independent parameter a 
in the profile function of Eq. ^ is replaced by a position-dependent effective one a'^^{Z) 
which evolves in a linear fashion along the formation length If from a reduced value for the 
SSC to the standard one associated with the normal hadron 

'Z < n^k? > [ ' 



a 



tot 
pN 



e{if- z) + e{z-if) 



(7) 



n \ If 

Here, n is the number of elementary fields (2 for the p"), kt = 0.350 GeV/c is the average 
transverse momentum of a quark inside a hadron, Z is the distance from the hard interaction 
point fi along the ejected hadron path, and Ti = is the hard-scale parameter that 
governs the CT effect. For the formation length If ^ 2p/AM^, we adopt the value AM^ = 



0.7 (GeV/ With this value we could reproduce the measured nuclear pion transparencies 



d. p decay The p° decays to a pair of pions with a branching ratio of almost 99%. For 
the kinematics of the JLab experiment the average life time of the p° corresponds to a mean 
path length of about 5 fm. This means that the majority of the will decay outside the 
nuclear medium and the anticipated effect on the computed nuclear transparencies is rather 
modest. The average opening angle of the pion pair is around 30 degrees in the lab frame. 
Therefore, it is fair to substitute the al% in Eq. g by {al"^^ + (7^+*^^) after the decay. The 
p° decay is expected to lower the nuclear transparency as the pion pair is subject to an 
increased attenuation compared to the p°. The adopted procedure gives an upper limit for 
the effect of p decay as it adds the attenuation on the two pions in an incoherent way [2] 



We include the effect of p° decay by replacing the pN total cross section a^°^ in Eq. Q by 
a position dependent one 



Here, Pp (Ep) are the rho meson momentum (energ y), and Tp = 149 MeV is the p° decay 
width in the lab frame. 

e. Cross section The cross section for p^ electroproduction on a nucleus A {e + A ^ 
e' + p^ + A*) takes the following form (Dirac spinors are normalized as uu = 1) 

^^""^ ^ \^ "EM Ee' A-f^'-^P (9) 

dE^,dVLe'dtd(t)p ^8(27r)4gEeQ2(i_e)^' ^' 

Here, the summation runs over all the shells a of the nucleus A, Ee and Eei are incoming and 
scattered electron energy, oem is the fine structure constant, the virtual photon q^{u, q) has 
four-momentum transfer = —q^ and degree of transverse polarization e, t = {q — PpY-, X] 
denotes the summing and averaging over spin degrees of freedom, and J^'^*^ is the matrix 
element of the 7* + A — > p + A* reaction (with A* the residual nucleus). In a factorized 
approach, the squared matrix element can be approximated as [211 [33] 

Y.\M^*^? ^ j rfp™^iAi^*^rpf'(Pr.) , (10) 

where MP'*^ is the amplitude of the 7* + — )■ p + process, pm = Pp — q is the missing 
momentum, and the distorted momentum distribution p^ is defined as 
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(11) 



In Eq.(lO), the matrix element can be related to the cross section in the center-of-mass 
frame: 

da'^*^ mlr v^..^*Ar f da 



d\t\d(t)* 87r2(s2 - 2s{m% - Q^) + {m% + Q2)2) . 

where s = W'^, the squared cm. energy for the 7*A^ system. In the last step we made 
use of the diffractive properties of the vector meson cross section at GeV energies. For the 
slope factor of the diffractive p° production, we take f3^p =6 (GeV/c)~2 [381 [i2] - [ii] and leave 



da 
d\t\ 



unspecified as it cancels in the nuclear transparency ratio. Combining Eqs. (9), (10) 



and (12) gives us the following formula for the cross section 
da''^ f ,^ OEM 



dEe'dQ^>dtd(f)p ^ J ^""'"4(271 



qE,Q\l - e)ml \d\t\ ' ' ^'"^ 







III. RESULTS 

The JLab Hall-B experiment E02-110 measured the nuclear transparency in p° electro- 
production for 0.8 < < 3 (GeV/c)2 on ^'^C and ^^Fe targets using the CLAS 47r detector. 
In order to select elastically produced mesons and suppress pions from resonance decays 
the following kinematical cuts were imposed: z{= Epjv) > 0.9 and W > 2 GeV. To ensure 
the selection of exclusive diffractive and incoherent events, t was limited to —0.4 < t < —0.1 
(GeV/c)2. In vector meson electroproduction, the virtual photon will fluctuate into a qq 
pair along a certain coherence length l^. = 2z//(Q2 -)_ m2-) and then scatter off a nucleon. 
To isolate a possible CT signal it is essential that Ic is more or less constant over all the 
kinematic ranges included in the analysis. As the qq is subject to initial-state interactions 
(ISI), a variation in Ic can cause a change in the transparency and thus mimic a CT signal. 
For the JLab experiment E02-110 one has 0.5 < /c < 0.8 fm. This is sufficiently smaller than 
the typical intra-nucleon distance in the nucleus, so ISI are not included in our calculations. 

The nuclear transparencies shown here are computed as the ratio of the cross section of 



Eq. (13) to its plane-wave approximation (PWA) equivalent. The PWA cross section neglects 



any form of FSI and is obtained by putting ^ = 1 in Eq. (11). It deserves highlighting 
that the experimentally extracted nuclear p transparencies use a ratio of the phase-space 
integrated cross sections from nucleus A to the deuteron. In Ref. [29j it is pointed out 
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that Fermi motion can introduce an overall change in the transparency with a constant 
factor and mimic part of a CT signal through a dependent change of the plane-wave 
ratio. The results of Ref. [22] suggest that in the PWA the ratio of the cross section on 
A relative to the deuteron increases with due to the applied cuts in t and the effect of 
Fermi motion. For the calculations presented here we prefer to compute the transparencies 
as a ratio of the RMSGA to the PWA cross sections for ^^C and ^^Fe. Indeed, the nuclear 
transparency computed as a RMSGA/PWA ratio for the same target nucleus is less prone to 
errors stemming from specific phase-space mismatches between ^^C, ^^Fe and the deuteron 
for example, as well as other uncertainties as possible medium-effects on the 7* + — > p + N 
reaction amplitude. 

The experimental kinematical cuts determine the phase space of the calculations. For 
each bin, we fix at its central value and compute the RMSGA and PWA cross sections 
for the ranges in (z/, z, t) used in the analysis of the data. The final result for the nuclear 
transparency at a particular value of the hard-scale parameter is computed as the ratio 
of the integrated RMSGA cross section to the integrated PWA one, whereby both in the 
nominator and the denominator the same ranges in {i>,z,t) are considered. 

Figure [2] shows the hard-scale dependence of the RMSGA nuclear ^^C(e, e'p") transparen- 
cies for cr*'^=25 mb, f3pN = 6 (GeV/c)~^, epN = —0.2 taken constant as a function of energy. 
We investigate the hard-scale dependence stemming from the effect of CT, p decay, and SRC. 
The computed transparencies that do not include CT show little dependence on Q'^ as can 
be anticipated from the energy independence of the Glauber scattering parameters at these 
energies. The SRC component in the FSI adds about 0.5% to the transparency, independent 
of the value of the hard scale parameter. Accordingly, the CT and SRC mechanisms can be 
separated by studying the hard-scale dependence of the transparency. The p decay lowers 
the transparency with about 7% at the lowest data point and 6% at the highest data 
point, reflecting the longer dilated p half-life time at higher energies. The CT effects induce 
the strongest dependence, increasing the transparency with 3 to 6% from low to high 
Q^. The inclusion of SRC in the FSI for the CT calculations yields almost no increase of 
the transparency as the produced p° in the SSC is already subjected to reduced interactions 
in the local neighborhood of the hard interaction where also the SRC effectively modify the 
density. 

Figures [3] and |4] compare the RMSGA transparencies including SRC and p decay with the 
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FIG. 2: (Color online) Relative effect of the SRC, of the p-decay, and of the CT on the predicted 
hard-scale dependence of the ^'^C{e, e' p^) nuclear transparencies. All calculations use (T*'jl^=25 mb 
and f3pN = 6 (GeV/c)~^. Base-line calculations (denoted as Trmsga) are obtained with the 
Eqs. ([l]) and ([2]). Dashed (sohd) curves include (exclude) the effect of p decay implemented with 
the Eq. ([s]). Black triangles (blue squares) exclude (include) the effect of SRC in the computation 
of the FSI. Green pentagrams include the effect of CT and red hexagons include both SRC and 
CT. Some points are offset on the x-axis for clarity. The kinematic cuts are those of the JLab 
E02-110 experiment (details in text). 

JLab data. Given the lack of detailed information on the cr*'^ parameter, we show results of 
calculations for a^°^=25 mb (Fig. |4]), and for o"*'^=20 mb (Fig. [3|. As mentioned in Sec. [ll] 
the calculations including the decay of the p meson represent a lower limit. We therefore 
have chosen to represent our calculations in Fig. |3]as a shaded region, confined on the lower 
side by the calculations with the p meson decay evaluated with the expression for the pN 
total cross section of Eq. ([s]) and on the upper side by the average of the calculations with 
p meson decay and those without. The width of the band reflects the estimated uncertainty 
in the computation of the effect of the p^ decay. In Fig. |4] we only show the curves of the 
calculations with the p meson decay in order to compare to alternate model calculations. 
The Figures [3] and |4] also include two other model calculations. The model by Frank- 
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FIG. 3: (Color online) Nuclear transparency for p electroproduction as a function of in ^^C 
and ^^Fe with a^°^ = 20 mb. Data are from [27j. Black curves are calculations from the model of 
Ref. [28j with (dashed) and without (full) CT effects. Green curves are calculations from the model 
of Ref. [29] with (dash-dotted) and without (dotted) CT effects. Blue (with CT effects) and red 
(no CT effects) shaded bands are results from the RMSGA model including rho decay and SRC 
effects. 

furt, Miller & Strikman (FMS) [28] is based on a semi-classical Glauber calculation and 
implements the effects of CT and p decay along similar lines as ours, with the values of 
(T*'J^=25 mb, AM2 = 0.7 (GeV/c^)^ and F^ = 149 MeV. The calculations of Ref. [29] from 
the Giessen group include a model for the elementary p production and describe the FSI 
with the semi-classical GiBUU transport model [lSllin|. Given their very different nature, 
it is satisfying that all three models yield similar ^^C transparencies when ignoring CT. The 
predictions including CT effects show a little more variation over the different groups, but 
display similar trends in their dependence. At identical parameter input in the Glauber 
part, we obtain a transparency that is about 5 % lower than the FMS model calculations. 
The difference between the calculations including and excluding CT effects is also bigger in 
the FMS model. It is worth noting that when considering the A-dependence, none of the 
models can satisfyingly describe the data for both nuclei with the same parameter set. This 
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FIG. 4: (Color online) As in Fig. |4]but for a^°^ = 25 mb. Blue circles (with CT effects) and red 
squares (no CT effects) curves are results from the RMSGA model including rho decay and SRC 
effects. 

was not the case for the A{e,e'p) and y4(e, eV"*") data, where one parameter set gave a very 
good agreement over the whole measured A-range. Comparing our calculations to the data, 
the RMSGA results for cr*'}!^=20 mb are a better match than those with cr*']!^=25 mb. 



In Table III we make a comparison between the computed and measured slopes of the 
hard-scale dependence of the transparencies. We clearly see a better agreement between the 
experimental slopes and those calculations including CT. As becomes clear from Fig. [3] the 
FMS and RMSGA calculations tend to underestimate the measured dependence of the 
nuclear transparencies. The GiBUU predictions have a stronger dependence. This can 
be explained by the effect of Fermi motion in large nuclei mentioned earlier in this section, 
which induces an additional rise in T with on top of the one due to the CT effect. The 
values of the parameters used for the p meson in the quantum diffusion CT model may be 
too conservative. Indeed, a bigger CT effect can be implemented with a smaller value of 
AM^ than the 0.7 (GeV/c^)^ adopted here. To date, however, there is little guidance with 
regard to realistic ranges for those parameters. Any further fine tuning of the parameters 
related to the magnitude of the CT effect is not opportune given the mentioned uncertainties 
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in the calculations. 



nucleus 


tot 

[mb] 


JLab data [27j 


RMSGA (+SRC 
+CT+decay) 


RMSGA 
(+SRC+CT) 


RMSGA 

(+SRC+decay) 


RMSGA 
(+SRC) 


12c 


25 


0.044± 0.015 ± 0.019 


0.029 


0.024 


0.0017 


-0.0038 


12c 


20 


0.044± 0.015 ± 0.019 


0.028 


0.021 


0.0032 


-0.0037 




25 


0.053 ± 0.008± 0.013 


0.020 


0.017 


0.0011 


-0.0050 


56pg 


20 


0.053 ± O.OOSib 0.013 


0.021 


0.016 


0.0026 


-0.0050 



TABLE I: Slopes [(GeV/c) ^] for the dependence of the different RMSGA calculations compared 
to the JLab data. 



IV. CONCLUSION 

We have performed relativistic Glauber calculations for A{e, e'p^) transparencies and 
compared them to the recent JLab data. The effects of short-range correlations, color 
transparency, and the decay of the p meson to a pair of pions can be included. The covered 
phase space in the calculations matches the experimental conditions and all kinematical 
cuts are taken into account. The predicted effect of SRC in the final-state interactions on 
the computed transparencies is small, or even negligible after also including CT effects. 
Including the p meson decay lowers the transparency up to 6-7 %, with a smaller decrease 
for higher Q^. This reflects the higher absorption rate of the two pions compared to the p. 

When comparing to the data and other model calculations, the results including the 
CT effect are consistently in better agreement than those without. The data suggest a 
stronger hard-scale dependence of the transparencies than predicted by our calculations 
using educated estimates for the parameters determining the magnitude of the CT effect. 
The presented comparison between calculations and data, gives additional support for the 
onset of CT in meson electroproduction reactions at energies of a few (GeV/c)^. 
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